A Pseudogap in the Single-Particle Density of States of a Tomonaga-Luttinger Liquid 
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We study a single-particle density of states (DOS) in a model of Tomonaga-Luttinger liquid with 
small (large) values of the Luttinger parameter Kc < 0.17, {Kc > 5.82), in the charge sector. We 
explain that such values of may be achieved by electron-phonon interactions without generating 
the spin gap. We suggest that electron-phonon interactions can be partially incorporated into the 
Tomonaga-Luttinger liquid scheme by introducing frequency dependent Kc{u!). We demonstrate 
that when the low-frequency asymptotic value -K'c(O) < 0.17 or K^{0) > 5.82, the single-particle 
DOS has a pseudogap behavior in frequency where p{uj)/uj vanishes at small frequencies. The DOS 
exhibits a peak the position of which scales as [-ftr(O) -|- {0)]ljo/2 where uo is the characteristic 
phonon frequency. 
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The Tomonaga-Luttinger liquid theory predicts a 
power-law frequency dependence of the single-particle 
density of states (DOS). In an apparent contradiction 
with this prediction many quasi-one-dimensional mate- 
rials (see, for example, ||l|) exhibit a pseudo-gap type 
structure of DOS. 

In this brief report we argue that it is perfectly possible 
to describe the pseudo-gap behavior in the Tomonaga- 
Luttinger framework provided one treats carefully the 
short-time part of the electron Green's function. 

Let us recall how DOS is calculated. The stan- 
dard non-perturbative technique used to study low- 
dimensional systems - the bosonization approach yields 
thermodynamic Green's functions in the Matsubara time 
representation G(r). To extract DOS one has to Fourier 
transform this function to get G(iw„) and then do the 
analytic continuation iw„ ^ to + iO. 

As we have said, in the limit T 0, the single-electron 
Green's function in the Tomonaga-Luttinger theory has 
a power-law asymptotics: 
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where A is a non-universal number and the exponent 
6 is determined by the interactions. For the standard 
spin-1/2 SU(2)-invariant Tomonaga-Luttinger liquid the 
exponent 6 is related to the Luttinger parameter Kc Q: 
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Assuming that G{—t) ~ —G{t) we can write down the 
Fourier transform as 
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dr sin(ti;„T)G(r) 
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At small frequencies this integral is dominated by the 
asymptotics (l|) when 9 < 1. Then we have 



G(z^„) « 2^Aj^ dr^-^^^ ^ (^„)' 
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The analytic continuation is quite straightforward in this 
case and one gets 
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When 6 > 1 the integral (||) converges at large Mat- 
subara times. However, at < 2 its second derivative 
with respect to a;„ is still determined by the asymptotics 
at large r such that we get 



G(icj„) = -iujng{0) - 2{uJn) 
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After the analytic continuation only the second term con- 
tributes to the imaginary part and we find that Eq. (|^) 
is still valid. Thus the value 9 — \ corresponding either 
to Kc = 3-2%/2 « 0.17 or K^ = 3-^2^2 w 5.83 marks a 
crossover into a pseudogap phase where p{uj)/uj vanishes 
at w ^ 0. 

Small values of Kc can be achieved in systems with 
strong retardation effects such as systems with electron- 
phonon interactions or Kondo lattices j^. In all these 
systems small values of Kc are achieved asymptotically 
at small frequencies. Therefore it is interesting to learn 
what is the area of validity of the universal power law 
behavior (||). Below we calculate the DOS for a model 
with electron-phonon interactions. 

The influence of phonons on the electron subsystem in 
quasi-one-dimensional metals have been studied by dif- 
ferent authors; mostly for the case of noninteracting elec- 
trons (see 0], 1^ and references therein). The combined 
effects of the Coulomb and electron- phonon interactions 
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was studied in J^, using the renormalization group 
approach. In j£| these effects have been studied in the 
framework of Tomonaga-Luttinger theory. 

Here we briefly repeat the derivation given in ^ . The 
lattice effects can be included in the Hubbard Hamilto- 
nian by making the hopping integral t dependent on the 
intersite distance: 

1 / 

(It has been proposed for the first time by Su, Schriefer 
and Heeger to describe the essential physics of conducting 
polymers |6)). Then the Hamiltonian takes the form: 
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where mj is dimensionless and k, has dimensions of energy. 
The Cj^a operators are the usual creation and annihilation 
operators for the electrons with spin a in the Wannier 
orbitals at site j and rij o- is the number of electrons. U 
is the repulsion of two electrons on the same site. 

In the case of an incommensurate band filling {Ak^ ^ 
27r/a) in the continuous approximation the electron- 
phonon part of this Hamiltonian generates a coupling 
between the lattice deformations and the 2k-p and Ak-p 
components of the charge density. The electron-phonon 
interaction contributes to an effectively retarded interac- 
tion between the electronic densities: 

S'int = - / drdr'dx ^ p{2lkF,x)Di{T - T')p{2lkF,x) , 

^ ] 1 o 
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where Di{t) is the phonons Green's function at q ~ 2lkp 
^ . Thus the main contribution to the interaction comes 
from phonons with large frequency. This interaction ef- 
fects both the spin and the charge sector. The phonons 
give a positive contribution to the current-current cou- 
pling constant in the spin sector (i.e. gs)- In this 
sector the electron-phonon interaction competes with re- 
pulsive forces responsible for go. We assume that the 
renormalized coupling constant is still repulsive such that 
there is no spin gap. 

In the charge sector the phonons influence the dynam- 
ics as well as the scaling dimensions, renormalizing the 
charge velocity Vc and Kc- In the limit of small frequen- 
cies ^ uji their renormalized values are: 
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where m* is interpreted as renormalized electrons mass 
with m being the bare mass. 



Since we are interested not only in the long time 
asymptotics of the Green's function, but in its inter- 
mediate time behavior, we shall model these behavior 
adopting the modified Gaussian model with the time- 
dependent Luttinger parameter: 
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where the function f{uj) takes values between /(O) = 
m*/m = and /(oo) = 1. In this brief report we 

suggest a semi-phenomenological form for the function 



fiiu) = 1 + 
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The above model yields the following Matsubara time 
single-electron Green's function 



G(r) 
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where ep is the ultraviolet cut-off (recall that we con- 
sider the Green's function at coinciding spatial points). 
Pictures of the single particle density of states are repre- 
sented on Figure 1. They are obtained by analytic con- 
tinuation of G{iu!n) from the imaginary axis to just above 
the real axis. 

The figures clearly show the crossover from the 
Luttinger- liquid type behavior with singular dp/duj at 
K > 0.17 to the pseudogap behavior at K < 0.17. An- 
other remarkable feature of these DOS is the peak at 
u! « K~^ujq/2 (see Fig. 2). Since the DOS is invariant 
under K this empirical formula can be general- 

ized as 
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It is interesting that the peak always occurs at frequen- 
cies larger than the characteristic phonon frequency uq. 
At small (large) values of K this discrepancy can be quite 
substantial. For example, a peak in DOS has been ob- 
served in K0.3M0O3 at w « 300 meV. The behavior of 
DOS at smaller frequencies is almost linear in lu which 
suggests K « 0.15. Then Eq. ( |l2| ) gives a reasonable 
estimate for the phonon frequency: Wq ~ 90 meV. 

E.P. is grateful to Steve Allen and particularly to Dave 
Allen for helpful discussions. 
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FIG. 1. Plots of the dependence of the density of states on frequency for different values of the parameters K and wo- 
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FIG. 2. Plots of the dependence of the density of states on frequency for different values of wo and fixed values of K. Here 
one can observe the placing of the peak of p(w) as function of wo. 
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